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ABSTRACT 

The physics of diffuse gas accretion and the properties of the cold and hot modes of 
accretion onto proto-galaxies between z = 2 and z = 5.4 is investigated using the large 
cosmological simulation performed with the RAMSES code on the MareNostrum su- 
percomputing facility. Galactic winds, chemical enrichment, UV background heating 
and radiative cooling are taken into account in this very high resolution simulation. 
Using accretion-weighted temperature histograms, we have perfomed two different mea- 
surements of the thermal state of the gas accreted towards the central galaxy. 
The first measurement, performed using accretion-weighted histograms on a spherical 
surface of radius 0.2 i? v ; r centred on the densest gas structure near the halo centre of 
mass, is a good indicator of the presence of an accretion shock in the vicinity of the 
galactic disc. We define the hot shock mass, Mshock, as the typical halo mass separat- 
ing cold dominated from hot dominated accretion in the vicinity of the galaxy. 
The second measurement is performed by radially averaging histograms between 
0.2i? v i r and i? V ir, in order to detect radially extended structures such as gas fila- 
ments: this is a good proxy for detecting cold streams feeding the central galaxy. We 
define M str eam as the transition mass separating cold dominated from hot dominated 
accretion in the outer halo, marking the disappearance of these cold streams. 
We find a hot shock transition mass of M s h oc k = 10 116 Mq (dark matter), with 
no significant evolution with rcdshift. Conversely, we find that Mstream increases 
sharply with z. Ou r meas urem ents are in agreement wit h the analytical predictions of 
Birn boim fc Dekel (|2003f ) and iDekel fc Birnboiml (|2006f ) , if we correct their model by 
assuming low metallicity 10~ 3 Zq) for the filaments, correspondingly to our mea- 
surements. Metal enrichment of the intergalactic medium is therefore a key ingredient 
in determining the transition mass from cold to hot dominated diffuse gas accretion. 
We find that the diffuse cold gas supply at the inner halo stops at z = 2 for objects 
with stellar masses of ab out 10 111 M(T ) , which is close to the quenching mass deter- 
mined observationally by iBundv et al.l (|2006f h However, its evolution with z is not 
well constrained, making it difficult to rule out or confirm the need for an additional 
feedback process such as AGN. 

Key words: methods: Numerical simulations, N-body, hydrodynamical, adaptive 
mesh refinement, galaxies: formation 



1 INTRODUCTION 

It is currently accepted that the ACDM theory provides 
a framework with which a large number of observed 
galaxy properties can be interpreted. This framework is 
referred to as the "hierarchical scenario of galaxy forma- 
tion" . Most importantly, this framework explains why many 
of these properties (physical sizes, black hole mass, bulge 



mass...) are found to correla te simply with galaxy mass 
(Kauffmann & Hachnclt 2000). Amidst this apparently sim- 
ple scaling of galaxy properties with mass, the discovery 
of a bimodality in the colour dis tribution of Sloan Dig - 
ital Sky Survey (SDSS) galaxies (|Kauffmann et al l 120031 ) 
stood unexpected and at odds with the predictions of hi- 
erarchical galaxy formation. Galaxy bimodality can indeed 
appear to be anti-hierarchical, as can be grasped from the 
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following simple argument. Under the assumption (com- 
mon in early semi-anlytical models, hereinafter SAMs) that 
the star formation rate (hereafter SFR) is proportional to 
the gas accretion rate, the lat ter being proportion al to 
the halo mass to some power (|van den Boschl [2002). one 
expects that, at all times, objects with the highest SFR 
should be the largest galaxies. In this framework, massive 
elliptical galaxies would still be blue and forming stars at 
z = 0. The fact that observed elliptical galaxies do not 
obey this fundamental prediction of the hierarchical sce- 
nario is the origin of the so-call ed "anti-hierarchical" be - 
haviour of massive red galaxies l|Rasera fc Tevssierl 120061) . 
This observation is further supported by the analysis of 
spectroscopic d ata, using star fo r matio n history reconstruc- 
tion methods ilReichardt et al.l [20011: iPanter et all 120031 ; 



ICid Fernandes et al.1 12004 lOcvirk et all l2006bl lal). Since 
these giant galaxies are in the form of apparently "dead" 
(i.e. no ongoing star formation) red elliptical galaxies, the 
quest has been ongoing for several years to find the origin 
of this halt in the star formation process (also refered to as 
"star formation quenching"). 

A substantial part of astrophysical research nowadays 
is devoted to searching for new physical mechanisms able 
to prevent cold gas accreted at i? v ir from falling into the 
galactic disc, condensing into molecular clouds and forming 
stars. The heating of the infalling gas from virialization and 
feedback from supernovae and hot stars has been considered 
as a serious candidate for more than a decade but seems 
insufficient to explain the drop in star formation o f massive 
systems in recent times (R asera fc Tevssierl 120061 ') , leading 
the authors to suggest a superwind phase for the high mass 
end of the galaxy population. 

Active Galactic Nucl ei (AGN) feedback has been pro- 
posed by several authors (|Bower et al.l 2006; Hopkin s et al.l 
120071 ) as the origin of this quenching, and has the additional 
desirable propery of preventing cool ing flows in the core 
of the most massive clust er galaxies l|De Lucia et al.l 120061 : 
ICattaneo fc Tevssierll2007l) . However, AGN physics are still 
poorly understood, both theoretically (because of the intrin- 
sic complexity of relativistic magn etohyd r odyn amic flows 
around black holes, see for instance iProgal (|2007l)). and ob- 
servationally (because of the small physical extent of the 
region of interest). Moreover, the mechanism involved in 
transferring the energy from the black hole accretion flow to 
the surroundings remains elusive. Jets have been proposed, 
and have the advantage of being supported by observations, 
but shock waves arising from the interaction of the jet with 
the interstellar medium would tend to push away the hot 
gas while leaving surround ing clumps of cold gas rather un- 
changed jSlvz et all 12005). However, this would depend on 
the position of the clumps with respect to the jet origin and 
the violence of the shock, and for instan ce, with a Mach 
numb er ~ 10 and density contrast ~ 10, [Nakamura et al.l 

(2006) do indeed predict cloud destruction. 

A jet-driven turbulence is another alternative, and to 
assess its relevance one has to repe at experiments such 
as tho se of iBaneriee et all (|2007l ) and lCattaneo fc Tevssierl 

(2007) at the galactic scale. The balance between the me- 
chanical power output of bubble-carving jets (estimated 
from radio luminosity at 1.4 GHz) and the radiative losses 
of the hot gas halo of galaxy clusters has been proposed 
as a signature of the global control of gas cooling by in- 



teraction with a central black hole jet (|Best et alj|2006l ). 
However, a necessary requirement for this scenario to work 
for galaxies is that the energy available in the jet remains 
in the galaxy, while observations show that radio sources 
have linear sizes significantly larger than their host galax- 
ies. It is thus not clear how jets can prevent star forma- 
tion, and it has also been argued th at they might actu- 
ally enhance star formation (|Silkll2005l ). More generally, the 
causal relation between AGNs and star formation is unclear: 
does the starburst trigger AGN activity and a subsequent 
quenching or does the AGN activity trigger the starburst? 
As a matter of fact, strong AGN activity is seen in galax- 
ies with intense star form ation activity (|Wild et al.l 120071 : 
Icid Fernandes et al.l l200ll ) , demonstrating that AGN and 
star formati on co-exist, although th i s could be just a short- 
lived phase dKauffmann et al. I l2003l : ISchawinski et al.ll2007l : 
ICiotti fc Ostrikerl 120071 ). Finally, purely radiative feedback 
from the ac cretion disc around the b lack hole might be an 
alternative (|Fabian et al.l [20061 . 120081 ). but the mechanical 
coupling through which the dust phase being blown away 
drags the cold gas along is uncertain. Hence, it is worth in- 
vestigating possible quenching mechanisms other than AGN 
feedback. In this respect, the detailed analysis of diffuse gas 
accretion around star forming galaxies is of great interest 
because it ca n provide a form o f self-r egulation. The sem- 
inal paper of iBirnboim fc Dekell l|2003l ) (hereinafter BD03) 
investigates the stability of hot accretion shocks around disc 
galaxies, showing that such shocks can exist only for haloes 
more massive than ~ 1O 115 M0. In an ideal spherical flow, 
this hot shock would prevent cold gas from reaching the 
disc (or at least sl ow it down) and t h us is likely to af- 
fect star formation. iDekel fc Birnboiml ([20061 ) (hereinafter 
DB06) extended this approach to the study of the stabil- 
ity of cold streams ("filaments") within the shock-heated 
halo gas. They showed that the observed transition mass 
from blue to red galaxies at 2 ~ could be matched to the 
critical mass at which a stable accretion shock can exist and 
that stable filaments would disappear around z — 1.5. These 
findings were also driven and further confirmed by numer- 
ical simulations of high redshift galaxy formati on based on 
smoot hed particle hydrodynamics (SPH), as in lKeres et al.l 
(2005). DB06 actually presented the rise of stable hot shocks 
not as the origin of the quenching but only as a necessary 
condition for an efficient AGN feedback. 

However, it is too early to discard the existence of sta- 
ble hot shocks or the destruction of the gas filaments as the 
origin of the galaxy bimodality. Indeed, no numerical study 
has considered the influence of chemical enrichment, which 
was shown in DB03 and DB06 to have a crucial impact on 
shock stability, since met allicity, along with gas dens ity, de- 
termines tlM_c£c4hig_mte jlSutherland fc Dopita|[l993l ). More 
recently, ICattaneo et alj (|2007 ) checked the good behaviour 
of SAMs with respect to hot and cold a ccretion modes 
by comparing Ga llCS (|Hatton et al.1 [2003) results to the 
iKeres et al.l l|2005h simulation, which does not take into ac- 
count chemical enrichment. 

In this paper, we address these very issues using the 
results of a very large cosmological simulation performed 
on the MareNostrum supercomputer at the Barcelona Su- 
percomputer Centre. It was performed within the Horizon 
collaboration (http://www.projet-horizon.fr) using the 
RAMSES code (Tevssicr 2002), which now includes a de- 
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tailed treatment of metal-dependent gas cooling, UV heat- 
ing, star formation, supernovae feedback and metal enrich- 
ment. 

The simulation parameters (Lbox = 50 /i _1 Mpc, 1024 3 
dark matter particles and a spatial resolution of about 1 
/i _1 kpc) are optimal to capture the most important proper- 
ties of gas accretion around typical Milky Way-like galaxies. 
The large box size allows us to have a large (more than 
100) sample of L* galaxies at z ^ 2 and above, with a 
strong statistical significance. For this mass scale, the main 
transitions in the accretion regime (appearance of hot ac- 
cretion shocks, disappearance of cold filaments) effectively 
takes place between 2 ^ z ^ 6. Moreover, the adaptive mesh 
refinement (AMR) method adopted in RAMSES allows us 
to investigate the flow in an Eulerian fra mework, in contras t 
to the Lagrangian approach adopted bv lKeres et ail (|2005l V 
Finally, it was recently observed that the appa rent bimodal- 
ity mass scale might increase with redshift l)juneau et al.l 
l2005l ; lBundv et aljbood : iHopkins et~aLll2007n , an intriguing 
behaviour that has yet to be checked against models. 

The outline of this paper is as follows: first we describe 
in Sec. 2 our methodology, in terms of numerical techniques 
and statistical measurements. We specifically introduce a 
new estimator to analyse the thermodynamical properties 
of accretion, namely accretion-weighted histograms. We then 
present in Sec. 3 our main results concerning the physical 
properties of the accreted gas. Our findings are then dis- 
cussed in the framework of earlier theoretical modelling in 
Sec. 4, and recent observations in Sec. 5. 



2 METHODOLOGY 

In this section, we first describe the MareNostrum simula- 
tion: a cosmological N body and hydrodynamics simulation 
of unprecedented scale with most of the physical processes 
involved in galaxy formation theory. We then describe our 
dark matter halo catalogue and its corresponding proper- 
ties, present our statistical tool - accretion-weighted tem- 
perature histograms, and discuss our criterion for separating 
diffuse gas accretion from satellite merging. 



2.1 The Horizon-MareNostrum simulation 

We have performed a cosmological simulation of unprece- 
dented scale, using 2048 processors of the MareNostrum 
computer installed at the Barcelona Supercomputing Centre 
in Sp ain. We have used intensively the AMR code RAM- 
SES l|Tevssieri I2OO2T ) for 4 weeks dispatched over one full 
year. This effort is part of a consortium between the Hori- 
zon project in France (http://www.projet-horizon.fr) 
and the MareNostrum galaxy formation project in Spain 
(http://astro.ft.uam.es/~marenostrum). This simula- 
tion should therefore be named the Horizon-MareNostrum 
simulation to avoid c onfusion with the GAD GET- 2 
MareNostrum simulation IGottlober fc Yepesl (|2007l ). In the 
rest of the paper, we refer exclusively to the Horizon- 
MareNostrum simulation unless explicitly stated otherwise. 

The main asset of this project relies on using a quasi ex- 
haustive number of physical ingredients that are part of the 
current theory of galaxy formation, and at the same time 
covering a large enough volume to provide a fair sample 



of the universe, especially at redshifts above one. Specifi- 
cally, we have considered metal-dependent cooling and UV 
heating using the Hardt and Madau background model. We 
have incorporated a simple model of supernovae feedback 
and metal enrichment usi ng the implementation described in 
Dub 01s fc Tevssierl(|2008h . For high-density regions, we have 
considered a polytropic equation of state with a 5/3 index 
to model th e complex, multi-phase and turbulent structure 
of the ISM (|Yepes et al.lll997l; ISpringel fc Hernquist|l2003r) 
in a simplified form (se e ISchave fe Dalla Vecchial 1 20071 ) ; 
iDubois fc Tevssierl i^POSt) ): the ISM is defined as gas with 
a density greater than no ~ 0.1 H/cm 3 . Star formation has 
also been included, for ISM gas only (nu > no), by spawn- 
ing star particles at a rate consistent with the Kennicutt 
law derived from local observations of star forming galaxies. 
Technically, we have p, — p gas /t, where t* = ^(nu/no)^ 1 ^ 2 
and to = 8 Gyr. Recast in units of the local free-fall time, 
this corresponds to a star formation efficiency of 5%. The 
simulation was started with a base grid of 1024 3 cells and 
the same number of dark matter particles, and the grid was 
progressively refined, on a cell-by-cell basis, when the local 
number of particles exceeded 10. A similar citerion was used 
for the gas, implementing what is called a Quasi-Lagrangian 
refinement strategy. Five additional levels of refinement were 
considered, but the maximum level of refinement was ad- 
justed so that the minimum cell size in physical units never 
exceeded one kpc. In this way, our spatial resolution is con- 
sistent with the angular resolution used to derived the Ken- 
nicutt law from observations. On the other hand, we are not 
in a position to resolve the scale height of thin cold discs so 
the detailed galactic dynamics are likely to be be affected 
by resolution effects. 

The simulation was ran for a ACDM universe with 
Q M = 0.3, Ov = 0.7, fl B = 0.045, H = 70 km/s/Mpc, 
as = 0.9 in a periodic box of 50 /i _1 Mpc. Our dark matter 
particle mass (mpart. ~ 8 x 10 6 M ), our spatial resolution 
(1 kpc physical) and our box size make this simulation ide- 
ally suited to study the formation of galaxies within dark 
matter haloes, from dwarf- to Milky Way-sized objects at 
high redshift. For large galaxies, we can nicely resolve the 
radial extent of the disc, but not its vertical extent, while for 
small galaxies, we can resolve the gravitational contraction 
of the cooling gas, but barely the final disc. 

A common goal of the Horizon collaboration and the 
MareNostrum galaxy formation simulation project is to in- 
vestigate the relative accuracy of the current SPH and AMR 
codes, through direct comparison of the results of a large 
galaxy formation simulation. To allow for this, the Horizon- 
MareNostrum simulation and the GADGET-2 MareNos- 
trum simulat ion use the same in it ial co nditio ns, which are 
descri bed in IGottlober fc Yepesl (120071 ) and iPrunet et~aH 
(2008). The simulation was stopped at redhift z ~ 1.5 be- 
cause the allocated time ran out. The total number of galax- 
ies at the end of the simulation was larger than 2 x 10 , the 
total number of star particles was more than 10 s , and the 
total number of AMR cells was larger than 5 x 10 9 . 



2.2 Virial spheres at rest 

In order to analyse the physical properties around high red- 
shift galaxies, we have built from our si mulation data a 
Friend-Of-Friend (FOF) halo catalogue (iEfstathiou et al.l 
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1 19881 ). For each snapshot and for each halo, we compute 
its mass and its centre of mass. The Virial radius is defined 
here as -R200&, the radius at which the average mass density 
in the halo is 200 times the background matter density. We 
then define the Bright Central Galaxy (BCG) of the halo as 
the highest gas density peak in a sphere centred on the halo 
centre of mass and of radius 0.57? v ir- This position will serve 
as our reference point for computing radial accretion rates. 
We checked a posteriori that this position also corresponds 
to the most massive substucture in the halo. We subtract 
from the gas velocity the mass-averaged velocity of the gas 
inside the Virial radius in order to put the system at rest. 
The basic physical properties of the gas are then mapped 
onto concentric shells centred on the BCG position. 



2.3 A new tool: accretion— weighted histograms 

It is quite common in cosmological simulations including gas 
physics to analyse the thermodynamical s tate of baryons 
using so-called "phase space diagrams" 1 l|Cen fc Ostrikei] 
ll9931 ; lKatz et al.lll996l ; iRasera fc Tevssierll2006l ). for which 
the total mass fraction in a given gas density and tempera- 
ture range is given as 2D histograms. This sort of diagram 
yields only a static view and does not include any reference 
to mass or energy fluxes. Instead, we propose to gain insight 
into the accretion regimes of cosmological haloes by using a 
new tool: accretion-weighted phase space diagrams. We still 
use temperature and density probability distribution func- 
tions, but we weight the contribution to each temperature 
and density bin by the local accretion rate. In this way, static 
regions will be discarded from the analysis, while large radial 
velocity regions will dominate the signal. Since the accretion 
is towards the BCG in the halo centre, we define concentric 
shells where the temperature, density, velocity and metal- 
licity are computed by smoothing the underlying 3D fields 
with a window function of scale R: 



T R (r, 



T(x')Wfl(x-x')d 3 x'. 



(1) 



We then sample this 3D field on various spherical surfaces of 
radius 0.2J? v ir < r < R v i T , with an angular resolution A9 ~ 
R/r. In the current analysis, we fix R ~ 2kpc, twice our 
spatial resolution, so that for our largest haloes, the sphere 
was sampled with 400 x 400 pixels. We will thus drop the 
subscript r for the rest of the paper. We obtained for each 
halo angular maps of each smoothed gas variables (density 
p, velocity v, temperature T, metallici ty Z) . We then defin e 
the local accretion rate m as (see e.g. lAubert et al.l (|2004h . 
lAubert fc Pichonl (|2007l )) 



m(r, Q.) 



dM 



— p v • nr 



(2) 



where the solid angle is defined by df2 = sin 6d9d(j) in the 
direction fl — (8, <f>), and n is the vector normal to the sphere 
and of norm 1. The total accretion rate across the sphere 
is recovered using M = J rhdfl. For most of this paper, we 
will omit the density variable, p, to focus on the thermal and 
chemical properties of the accretion flow, described here by 
T and Z. At a given radius, r, we can marginalize Eq. (J5J) 



1 as in different physical phases, not position- velocity phase space 



over all cells which have a given temperature, T, and obtain 
the accretion rate per unit temperature as 



f dM 
m(r,T) = J S D (T-T(n))mdn= — . 



(3) 



where Sr> is the Dirac function. The total accretion rate 
across the sphere is recovered using M = J rhdT. Simi- 
larly, while marginalizing over all angles which have a given 
temperature, T, and a given metallicity, Z, we may intro- 
duce our main tool, the accretion-weighted temperature- 
metallicity two-dimensional probability density function 
(hereinafter PDF) 



(r,T,Z)= S D (T - T(n))5 D (Z - Z{Q,))mdQ. = 



d 2 M 
dTdZ 



(4) 



2.4 Diffuse accretion versus clumpy satellites 



In this paper, we are interested in characterizing the ac- 
cretion of diffuse intergalactic gas rather than the accre- 
tion of galaxy satellites. We therefore need to separate the 
contribution of infalling gaseous discs from the smooth ac- 
cretion through filaments or other diffuse components. In 
our model, the star-forming dense ISM is defined as nu ^ 
0.1H/cm 3 . We remove from our spherical analysis all pixels 
whose density exceeds this theshold. Although this trun- 
cation may seem brutal, infalling satellites are frequently 
embedded in diffuse filaments. As such, the filament is the 
natural surrounding of the infalling satellite, and deciding 
where the boundary lies is difficult yet critical. In this con- 
text, a density criterion is still the most straightforward sep- 
aration method. In the future, a possible alternative could 
be to use a segmentation algorithm to separate the filament 
from the satellites on topologi cal grounds, in the spirit of 
the skeleton reconstructions of ISousbie et alJ (|2008l ). 



2.5 Entire halo versus galaxy vicinity 

We use two distinct estimators, m and (m), to characterize 
the diffuse gas accretion around the central galaxy. 
The first one, rh is directly related to the gas properties 
close to the galactic disc: we compute the accretion-weighted 
T — Z histogram, Eq. Q, at radius r = 0.2 i? v ir, which 
turned out to be close enough to, but not intersecting, the 
neutral HI disc. This region defines the disc vicinity, for 
which we would like to analyse the thermal properties of the 
accreted gas. We easily detect the presence of an accretion 
shock upwind of the disc if the accreted gas is predominantly 
in a hot phase. On the other hand, if the accreted gas is pre- 
dominantly in a cold phase, it means that no accretion shock 
is present above 0.2i? v i r . 

The second estimator, (rh), is based on averaging the 
accretion-weighted histograms measured at different radii 
between 0.2i? v i r and i? v ir-l 



(m)(T,Z) 



dv 

m{r,T, Z) — . 

Ar 



(5) 



where Ar = 0.8 R v n and rh(r,T,Z) is given by Eq. Q. 
Coadding different histograms at different radii increases 
the weight of coherent radial structures such as filaments 
or cold streams that eventually extend out to (or beyond) 
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the Virial radius. Indeed, if the velocity flow is chemo- 
thermodynamically similar between two neighbouring shells, 
the corresponding histograms will add up consistently. Note 
that under the assumption of a steady spherically symmetric 
flow, the accretion rate does not depend on radius. As we will 
see later on, even in the presence of an accretion shock, cold 
streams can persist in the halo and directly feed the central 
galaxy with fresh cold gas. This phenomenon will appear in 
our histograms as a dominant cold phase. If, on the other 
hand, these filaments are destroyed, the histograms will be 
dominated by the hot phase out to the halo Virial radius. 



3 PROPERTIES OF DIFFUSE GAS 
ACCRETION 

We computed the accretion-weighted PDFs for several hun- 
dred haloes spanning dark matter masses between 10 10 — 
10 13 M Q between 2 < z < 5. We then co-added (stacked) 
these distributions for haloes of the same mass range in or- 
der to produce an "average" PDF for a given mass scale. 
We use these stacked PDFs to study the typical tempera- 
ture and metallicity distribution of the accretion flow, and 
the transition mass between the hot dominated and the cold 
dominated accretion regimes. Throughout the whole paper, 
we use exclusively a basis 10 logarithm, log 10 , and will there- 
fore drop the subscript io in text, equations and figures. 



3.1 Bimodality in the temperature distribution 

3.1.1 Hot and cold modes 

The left and middle columns of Fig.[T]show several radially- 
averaged accretion-weighted stacked-histograms (following 
Eq. ©) for haloes from 10 10 to 10 13 M Q taken from the z = 4 
and z = 2.5 snapshots of the simulation, respectively. The 
right panel shows these histograms for the 0.2_R v i r spherical 
surface at z = 2.5. The number of galaxies in each mass bin 
is given in the top right corner of each panel. These numbers 
reflect the halo sampling strategy and not the halo mass 
distribution of the simulation. For lisibility, the histograms 
have been convolved with a smoothing kernel of about 0.2 
dex FWHM in T and Z. We see that the accretion pattern 
involves two main distinct components: 

(i) A cold component, the metallicity of which extends 
over several decades 

(ii) A hot, relatively metal-poor component, the temper- 
ature and contribution of which increases sharply with halo 
mass 

As such, accretion itself is clearly bimodal in temperature. 
Indeed, at any halo mass, little mass is ever accreted around 
T ~ 1 — 2.5 x 10 s K. Instead, most of the mass is accreted 
either belo w or above t h is tem perature. As already noted by 
BD03 and [Keres et~ail l|2005h . this involves a link with the 
physics of cooling: little mass will be accreted at tempera- 
tures where the cooling is efficient, since gas cannot remain 
at this temperature for very long. This provides a natural 
temperature threshold that allows us to separate the cold 
and the hot accretion modes, associated with low and high 
mass haloes respectively. The middle and right columns of 



Fig. [T] show that this bimodality in temperature is also seen 
at z = 2.5, in the whole halo as well as at the inner halo. 

A good criterion for establishing the existence of a well- 
developed hot phase could be to require the existence of a 
saddle point in the accretion-weighted PDFs of Fig. [1] as 
a local minimum in temperature and a local maximum in 
metallicity. According to this criterion, the hot phase of the 
lowest mass bin haloes is not well-developed at any redshift 
between z = 2-4. The accretion-weighted PDFs can be nor- 
malized using the total diffuse gas accretion at the Virial 
radius. Fig. [2] shows that the accretion rate decreases with 
cosmic time for all masses and increases with mass at fixed 
redshift. This evo l ution is very similar to t h at reported in 
iKeres et all ll2005l):lRasera fc Tevssierl (|2006| ); IGuo fc White] 
1 20071 ) : lNeistein fc Dekell i|200Sl ). in trend and normalization. 



3.1.2 Identification of the phases 



The cold component described in Sect. [5TT7T1 can itself be 
decomposed into two distinct sub-components with differ- 
ent behaviours, based on their metallicity. On one hand, the 
high-Z part of the cold component is made of the surround- 
ings of satellite galaxies, i.e. a region close enough to the 
satellite core to be chemically enriched, but far enough to 
escape the density truncation described in Sect. 12.41 This 
can be clearly seen on Fig. [5] This component is seen at 
all redshifts. Conversely, the low-metallicity tail (Z ^ 10 -3 
Zq) of the cold component is very prominent at z = 4 but 
disappears at z = 2.5 in the inner and outer halo. It can be 
identified with the dense, cold, metal poor gas filaments seen 
in Fig. [S] This is best seen at z — 4, where the filaments are 
better defined. Their density is somehow intermediate be- 
tween the galaxy discs and the background. The left panel 
of Fig. [S] clearly shows a cold metal-poor filament tunneling 
down all the way from almost 2 i? v ir to the central galaxy's 
disc, while the cold metal-rich denser phase makes up the 
surroundings of satellites. The hot phase, with intermediate 
density and metallicity, is distributed in a large bubble, of 
radius smaller than R v i r at z — 4 but significantly larger 
than ii v i r at later times. 



3.2 Metallicity of the hot mode 

3.2.1 The heterogeneous hot phase 

The top panels of Fig. Q] show that the metallicity of the hot 
and cold phase can differ by up to 3 decades. This gap in- 
volves a huge difference in the ability of the gas to cool down 
radiatively, further aggravated by the relative densities of 
the two phases. This bimodal metallicity distribution is also 
seen at the inner halo. This has important consequences for 
models of galaxy formation and evolution, and highlights 
the necessity of treating the gas as being composed of two 
phases with widely different metallicities in SAMs. More- 
over, at a given temperature, the large spread of the PDF in 
metallicity shows that the gas is not well mixed. Indeed, a 
perfectly well mixed accreted gas would show up as a narrow 
peak in metallicity. This non- homogeneity is not a result of 
the stacking of the PDFs of haloes with different metallici- 
ties, but is already seen in individual haloes, although in a 
noisier fashion. Neither is it the result of stacking the PDFs 
of several shells of different radii: the top right panel shows 
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z=4 [0.2-1] FU z=2.5 [0.2-1] Rvir z=2.5 0.2 R vi , 




z = 2.5, radially averaged. Right: z = 2.5, 0.2 R v i r . The numbered labels on the contours give the logarithm of the PDF. The number of 
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Figure 2. Diffuse gas accretion rate measured at i? v ir versus dark 
matter halo mass at various epochs between redshift 2 and 5.4 



that the hot phase metallicity distribution is already 2-3 dex 
wide in a single shell at the inner halo for the most massive 
halos. 



3.2.2 Accretion-rate weighted metallicity 

Having acknowledged the chemical heterogeneity of the hot 
phase and its possible implications, we now turn to defining 
the mass accretion-rate weighted metallicity Zhot of the hot 
phase, as the metallicity the hot gas would have if it was 
perfectly well mixed. 



Zhot (r) 



f™rh(r, T, Z)ZdTdZ 
m{r, T, Z) AT&Z 



(6) 



Fig. [3] shows the variation of Zhot with respect to radius for 
several mass bins at z = 2, and for the largest mass bin at 

2 = 3. 



3.2.3 Metallicity profile 

For all mass masses, Zhot increases sharply towards the halo 
center, and all the haloes have a similar central metallicity 
Zhot(0.1 -Rvir) = 0.1 Zq. However, the slope of the metal- 
licity profile appears to depend strongly on mass, and the 
outskirts of the more massive haloes are found to be more 
metal-rich than their low mass counterparts. It is not clear at 
this stage of the analysis of the simulation if this trend can 
be explained by metal-rich winds from the central galaxy 
alone or stripping and dilution of the metal-rich disks of 
infalling satellites. It is also likely that winds in infalling 
satellites and stripping work together to enrich the halo 
gas: winds deposit metals out of the plane of the satellite's 
disk, and hence facilitate the stripping of these metals from 
the satellite. 



We also show that these metallicity profiles are well 
fitted between r/R vir = 0.1-1 and Mdm = 10 10 - 10 13 M Q 
by the simple logarithmic law: 



Zhot(r) = a + b log(r/i? vir 



(7) 



where the dependence with respect to mass is restricted to 
the coefficients a and b given by: 



a = -8.725 + 0.5 log(M DM ) 
b = -7.725 + 0.5 log(M DM ) 



(«) 
(9) 



The thick gray lines in Fig. [3] show this fit for Mdm = 
10 12 ' 5 M Q (upper line) and Mdm = 10 11 ' 5 M (lower line). 



3.2.4 Evolution with redshift 

Fig. [3] also shows the metallicity profile of the highest mass 
bin for z = 3. The latter profile seems to be a mere down- 
wards translation of the z — 2 high mass profile, indicat- 
ing that the enrichment rate s — —d\og(Zhot/ZQ)/dz w 
0.2 — 0.3 is almost the same at the inner halo and at the outer 
halo. This compares rather well wit h the average chemica l 
enrichment rate s — 0.17 found by |Pe Lucia et all l|2004t ). 
It is also in agreement with t he central chemical enrich- 
ment rate of ICora et all (j2008T ). who use a semi-analytical 
approach. 



3.2.5 X-ray clusters 

To roughly compare our profiles with metallicity measure- 
ments of X-ray clusters, we first need to extrapolate linearly 
the evolution of Zhot at the inner halo down to z = 0. 
This gives log(Z ho t(0.1i?vir, z = O)/Z ) m —0.4 (assum- 
ing s = 0.3), which is compatible with observations of 
IVikhlinin etahl (|2005l ). Indeed, they find log(Z/Z ) = -0.2 
at 0.1 R.500- The agreement further hol ds when comparing to 
the low-redshift part of the clusters of iBalestra et al.1 (|2007T l 
and iMaughan et al.l (|200&t ) . Note however that the galaxy 
clusters studied in the latter papers have progenitors more 
massive than the most massive halos of our simulation, and 
thus expectedly overall higher metallicity. Indeed, our most 
massive halos are more likely to end up in groups of galaxies 
at z — 0. Such objects also h ost hot X-ra y emitting gas with 
about half-solar metallicity (jBuotc 200Cf), which is in agree- 
ment with our extrapolation. Our finding that the metallic- 
ity we measured at O.liivir does not depend on halo mass is 
further supported by the fact that observed galaxy clusters 
and galaxy groups have roughly similar central metallicity 
(while the metallicity at the outskirts can be very different). 
Finally, the enrichm ent r ates found for X-ray c lusters by 
IBalestra et ail l|2007l ) and IMaughan et~ai1 (|2008h (s « 0.3 
from z = to z = 1.3 are well matched by the enrichment 
rate measured in the simulation between 2 = 2 — 3. 



3.3 Two critical masses for diffuse gas accretion 

Marginalizing the accretion rate over metallicity and inte- 
grating over temperature on the hot and cold temperature 
domains yields the hot and cold accretion rate respectively. 
Dividing by the total accretion rate at the chosen radius 
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Figure 3. Metallicity of the hot accreted gas as a function of 
the normalized radius r/R vlI at z = 2 for halo mass between 
10 10 - 1O 13 M , and z = 3 for the highest mass bin only. The 
thick gray lines shows the analytical fit of Eq. [7] for AfrjM = 



10 12 5 M (upper line) and Mdm 



10 1 



1 M (lower line). The 



dependence with respect to radius is very strong. While the profile 
gets steeper with decreasing mass, the innermost metallicity is 
almost constant. 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



hot and cold fractions 

xcmxnz 




l_ [0.2-1 ]R 

m i i i M^mrffTTiTT 



U.LLJJjL^ 




0.2 R, 



10.0 10.5 11.0 11.5 12.0 

l0g(M D M[M S un]) 



12.5 



Figure 4. Evolution of the hot (thin line) and cold (thick lines) 
accreted gas mass fractions versus Mdm for z 6 [2,5.4]. Top: 
(/cold) an d (/hot) (integrated from 0.1-f? v i r down to 0.2i? v j r ). Bot- 
tom: / co id and /hot on the 0.2i? v i r sphere. 



gives the contributions of the hot and cold mode to the to- 
tal accretion rate: 



/coid(r) = 



/hot(r) = 



T=T rZ = oo 



M(r) JT=0 JZ=0 
rZ = oo 



m(r,T,Z)ATAZ , (10) 



M(r) Jt=t Jz=o 



m(r,T,Z)ATAZ , (11) 



where rh(r,T,Z) is given by Eq. A similar expression 
involving (m(T,Z)) and Eq. <(Sj allows us to define (/cold) 
and (/hot)- The top panel of Fig. [3] displays the fractions 
computed from the accretion-weighted histograms averaged 
over the entire halo (between 0.27? v i r and R v i T ). The bot- 
tom panel of Fig. U shows these fractions measured at ra- 
dius 0.2i? v ir (inner halo) as a function of mass for various 
redshifts. A common feature of these plots is the increasing 
importance of the hot accretion mode with increasing mass, 
and the corresponding decreasing contribution of the cold 
mode, as could be foreseen from the top panels of Fig. [T] 
The mass at which (/ co id) = (/hot) defines the critical mass 
marking the transition between the two accretion regimes. 

This critical mass seems to increase sharply with red- 
shift (radially averaged case, top panel of Fig. Note that 
at redshift 5.4, it can only be guessed since no halo in the 
simulation is massive enough to have (/hot) *5 0.5. This evo- 
lution is the signature of a gradual disappearance of cold ra- 
dially extended features, like filaments, in the massive haloes 
between z = 5.4 and z = 2. This is illustrated by Fig. [5] 
showing maps of a typical halo of mass Mdm = 2 x 1O 12 M 
at z = 4 (left) and another halo of the same mass at red- 
shift z =2 (right). While the former features clear filaments 
streaming into the inner halo, the latter lies at the centre of 
a hot bubble, with no apparent filaments inside the Virial 
radius (large black circle). The critical mass defined by the 
accretion transition in the entire halo marks the disappear- 
ance of cold streams. It is therefore called here M a t re am. 

On the contrary, the critical mass at the inner halo 
(0.2i? v i r ), i.e. the mass of halos with /hot ^ 0.5 shows only a 
slow variation with redshift, if any. It indicates that accre- 
tion in the inner parts of the halo switches to the hot mode 
as soon as Mdm *5 1O 11,5 ~ 12 M at all redshifts considered, 
while the outer part of the halo can still be dominated by 
the cold mode. Again, this is well illustrated by Fig. [5] The 
inner part of the halo (inner circle) is shock heated at both 
redshifts, although the radius of the accretion shock is much 
larger in the low redhift snapshot. At high redshift, the ac- 
cretion shock coexists with cold streams coming from the 
outer parts of the halo. The critical mass defined at the in- 
ner halo marks the appearance of an accretion shock around 
the galaxy. It is therefore called here M s hock- 

An important parameter in our approach is the density 
threshold we used to remove clumpy satellites from the anal- 
ysis. We checked that our conclusions in terms of transition 
masses and average metallicity of the hot phase are robust 
to changes of this parameter by repeating our measurements 
with a lower density threshold. The main effect of this extra 
gas removal is to reduce the fraction of high metallicity, cold 
gas in the vicinity of galaxy satellites, without any notice- 
able effect on the metallicity of the hot phase and on the 
critical masses. 
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Figure 5. Maps of the physical properties of the gas for a typical M = 2 X 1O 12 M0 halo. Left: z = 4, right: z = 2.5. Top: Projected 
density, middle: temperature slice, bottom: metallicity slice. The large (small) black circle shows the location of iJ v j r (0.2H v i r ) for the 
central galaxy. Note the disappearance of filaments at low redshift. 
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4 COMPARISON TO EARLIER 
THEORETICAL MODELLING 

The physics of accretion has been investigated by several 
authors in the past. It is insightful to review their results 
in light of our measurements. Fig. [6] shows the evolution 
of Mstream with redshift. For z = 2 it is approximately 
10 11,5 Mq and it increases sharply with increasing z. At z ^ 4 
and above, even the most massive haloes in the simulation 
are still dominated by cold accretion. However, a rough ex- 
trapolation of the 2 = 4 curves of Fig. [4] yields a transition 
mass of about « 1O 12 ' 7 M0. Qualitatively, this behaviour is 
in agreement with the evolution of Af s t rC am with redshift as 
derived by DB06. According to their study, filaments exist 
only in M ^ M stre am haloes. However, their figure 7 shows 
that M s treamDB06 (z = 2) = 10 12 ' 4 Mq , which is significantly 
larger than what we find. We will now show that correct- 
ing the metallicity assumptions of DB06 with our measure- 
ments (low metallicity in the filaments) can reconcile these 
discrepant values. 



4.1 Approximation of DB06 model 

In order to compare our results with DB06 calculations, we 
will first build a crude interpolation of their M stre am and 
Afahock for any redshift, radius, and metallicity. Analysing 
the dependence of M s hock with respect to the problem pa- 
rameters in equation (34) of DB06, we can isolate the de- 
pendence in metallicity as: 



log(M sh ock(r,Zo,z)) = 0.71og(Z ) + A(r,z) . 



(12) 



where Zo = Z/Zq at z = 0, and A(r,z), which contains the 
dependence of M sriOC k with respect to radius and redshift, 
can be tabulated from fig. 2 and 4 of DB06. Then, taking 
the log of their equation (40), we see that M B tream is defined 
with respect to the critical mass for shock stability, M s hock, 
as 

log(M stroam ) = 2 log(M shock ) - log(M») - log(3) , (13) 

where M* is the typical dark matter halo ma ss at a given 
redshi f t as computed usin g the formalisms of Lahav et al.l 



[|l99ll ); ICarrorl et all (|l992l ); lMo fc White] l |2002j ) and follow- 
ing Appendix 3 of DB06. Since Zo is required at z = 0, we 
need to extrapolate our measurements at z — 2 to z = 0. To 
do so, we will for consistency use t he same enrichment ra te 
s=0.17 as in DB06, as derived from lDe Lucia et all l|2004h . 



4.2 Hot shocks 

Since the volume filling factor of the cold phase is in gen- 
eral negligible compared to that of the hot phase when 
looking at the halo as a whole, only the hot phase metal- 
licity Zhot is relevant when evaluating the stability of the 
hot shock. At the inner halo, Fig. [3] or the fit of Eq. [7] 
show that Zhot(0.1R v i T , z = 2) ~ — 1. This translates into 
log(Zhot(0.1-R v ir, z — 0)/Zq) — —0.65, which is rather close 
to the metallicity assumption of DB06, i.e. log(Zo/Z©) = 
— 1 at 0.1 i?vir- As a consequence, we expect the criti- 
cal hot shock mass M srioc k computed by DB06 to match 
our measurements well, provided that we compute and 
measure M s h oc k at the same depth in the halo, i.e. at 
0.2 R v ii. We take as an average metallicity at this radius 



log(Zhot(0.2J? v i r , z = 2)/Zq) = —1.7, which is intermediate 
between the 3 mass bins shown in Fig. [3] This translates 
into log(Zhot(0.2_R vir , z = 0)/Z e ) w —1.4. The shock mass 
Afshock obtained with this metallicity and radius is shown 
in Fig. [6] A good agreement is indeed achieved between 
this model and our measurements . We also agree w ith a 
quasi constant M ahock as found bv lKeres et alT f2005) and 
iBirnboim et al.1 (|2007h in SPH simulations, although the ab- 
solute normalization seems to differ. However, their method- 
ology is quite different (they analyse the temperature history 
of their gas particles) and their metallicity is uniform and 
constant throughout the whole simulation. 



4.3 Cold streams 

In DB06, Afstrcam is related to M sri0 ck, which strongly de- 
pends on metallicity via the cooling function, as also shown 
by fig. 10 of BD03 and fig. 2 of DB06. As established from 
Sect. l3.ll and Fig.[T] the metallicity of cold streams is rather 
low. At z = 4, for example, the metallicity queue of the cold 
accretion mode extends down to Z = 10~ 8 Zq and possibly 
lower, while the cold high-Z accretion consists of the sur- 
roundings satellite galaxies gas disks, rather than genuine 
cold streams. The efficiency of radia tive cooling decreases 
towar ds low metallicity, and fig. 13 of ISutherland fc Dopital 
(1993) shows that at Z/Zq — 10 -3 the cooling properties 
of the gas are already those of a primordial mixture. We 
use this value of the metallicity for the DB06 modelling of 
Mstream- Doing so shifts the disappearance of cold streams 
to earlier times (i.e. higher z) with respect to the original 
assumption of DB06 (Zo = O.IZq), as a consequence of less 
effective radiative cooling. Fig. [6] shows that this modelling 
of the low-Z cold streams then match our measurements 
rather well. The value plotted for M strcam at z — 5.4 is an 
extrapolation and is given only as a lower limit (indicated 
by the arrow). 



4.4 Discussion 

We see that while the relatively high metallicity of DB06 
at the inner halo is a fair assumption, a hundredfold lower 
metallicity must nonetheless be assumed for the gas fila- 
ments, which then have radiative properties similar to a 
primordial mixture. Once the importance and the effect 
of these assumptions has been accounted for, the consis- 
tency between the MareNostrum measurements and the the- 
ory of DB06 is remarkably good, and shows that their an- 
alytical approach indeed seems to capture the essence of 
our (nonetheless limited) understanding of the processes in- 
volved in gas accretion physics as modelled by the simula- 
tion. 

The agreement on M str eam also suggests that the main 
process driving the stability of cold filaments in a hot halo is 
the competition between a compressive increase of temper- 
ature in the filament on the one hand, and radiative cooling 
on the other hand. This is an important point since sev- 
eral other hydrodynamical processes are expected to be at 
work, such as Kelvin-Helmholtz instabilites, which are not 
modelled in DB06. However, th e MareNostrum si mulation 
is not tailored to resolve them. iKlein et al.l ()l994l ) showed 
that at least 100 cells per cloud radius are needed to resolve 
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Figure 6. Evolution of Af s hock and Afstrcam with redshift, from 
our measurements and comparison to analytical modelling. The 
solid line shows DB06 prediction for M s h O ck(0-2i? v j r ) assuming 
Zo/Zq = 0.04, while the dotted line shows their prediction for 
-^stream with a modified metallicity Zo/Zq = 10~ 3 . Finally, the 
dash dotted line shows the constant transition mass reported by 
Keres et al. (2005). The error bars on our measurements are equal 
to the widths of the mass bins used. The arrow pointing to the 
Mstream (z = 5.4) indicates that this point is given as a lower limit 
only, as can be estimated from Fig. 4. 



interface instabilities involved in cloud destruction by shock 
waves. However, it has been argued that the gas filaments 
can be stabilized by the underlying dark matter stream, and 
could totally prevent Kelvin-Helmholtz instability from ap- 
pearing. Filaments can also be stabilized if the flow is su- 
personic, which is indeed the case: for the halo shown in 
Fig. the inflowing filament from the lower right quadrant 
of the central galaxy has a Mach number of 2-3. Moreover, 
the density contrast of the cold gas stream and the average 
density at the Virial radius is also an important parame- 
ter. DB06 note that it is not a well-constrained quantity 
at the moment, although measurements are under way, and 
could well shift within a whole decade around the current 
assumptions in their modelling. Decreasing this ratio could 
have a similar effect to decreasing the metallicity, in making 
the cold filaments disappear earlier (i.e. denser filaments are 
more stable). 

Finally, we recall that the low Z measured in the fila- 
ments could be a result of the limited resolution of the simu- 
lation. Indeed, additional fragmentation down to the lowest 
masses allowed by the UV radiation field could still take 
place in the filaments, forming high-redshift dwarf galaxies, 
and the resulting star formation episodes could enrich the 
filaments. 



5 COMPARISON TO OBSERVATIONS 

It is difficult to find observables to compare our measure- 
ments to since our measurements "only" reach down to 
z = 2 while observational studies of the galaxy bimodal- 
ity are generally restricted to z ^ 1.5, and the bimodality 
becomes clear only at later times. Also, the definition of the 
critical masses constrained by theoreticians and observers 
have little in common: theoreticians have access to quanti- 
ties that are in general impossible to measure with existing 
or even future observing facilities (for instance the mass ac- 
cretion rate of cold gas at R v i T ), while observers have access 
to enormous volumes of the universe and events which would 
be impossible to simulate with the sufficiently high spatial 
and temporal resolution required for imposing useful con- 
straints. Hence, shortcuts are taken. They always involve 
external assumptions, and the relevance of the quantities 
being compared must be carefully examined. Here we dis- 
cuss our results in the light of t he observed t r ansiti on and 
quenching masses measured by iBundv et al.l (|2006h , here- 
inafter B06. 

The stellar to dark matter mass ratio M*/Mdm used to 
estimate the observed stellar mass M* from the dark mat- 
ter halo mass Mom and its evolution with redshift is sub- 
ject to large uncertainties and is a hotly debated topic. 
It is linked to the evolution of the Tully -Fisher relation 
l|Tullv fc Fisherlll977l : [Beil fc de JonefeOOlf ). However we re- 
call that the latter rather focuses on the M*/M tota i ra- 
tio where the total mass Mtotal = M + M gas + Mum 
is the sum of the stellar, gas and dark matter masses. 
However, the sum is always strongly dominated by the 
dark matter mass for the systems and radii we con- 
sider here. It seems reasonable to expect that the ratio 
M*/M DM increases with decreasin g redshift ( Rettura et al.l 
120061 : iKannappan fc Gawiserl 120071 '). as more gas is turned 
into stars. However, a non-evolving M * /Mum could also be 
in agreement with bo t h observations ( Bamford et al.l 120061 : 



iBohm fc Ziegler 20071: Atkinson et al. 20071 ) and numerical 
studies ( Portinari fc Sommer-Larsenl |200*7T ). In our simula- 
tion a ratio M* /Mum ~ 50 seems reasonable at z — 2 for 
structures in the mass range considered here. 
BOG measured the evolution with redshift of the transitional 
mass Mtr marking the transition from the blue to red se- 
quence of galaxies. This study is based on a sample of 8000 
DEEP2 galaxies between 0.4 ^ z ^ 1.4. The galaxies are 
divided into two groups (star forming and passive), ideally 
according to their SFR, with a limiting SFR=O.2M /yr. 
Howe ver, the [OH] emissi on line generally used to derive the 
SFR dKewlev et al.ll2004l ) falls in the DEEP 2 survey wave- 
length range only for z > 0.75. As a consequence, the colour 
index (U-B) is also used as a proxy for the SFR as measured 
from [Oil] , and an absolute magnitude-d ependent colour cut 
is ma de as an alternative to the SFR cut l|van Dokkum et al.l 
2000). The authors show that indeed, very few galaxies with 
(U-B)^ 0.2 are forming stars. However, there is a large frac- 
tion (« 30%) of galaxies with (U-B)$; 0.2 which are passive. 
This means that the sample of passive galaxies will suffer 
only limited pollution from star-forming galaxies, while the 
star-forming sample contains a significant fraction of passive 
galaxies (~ 30% at z — [0.75 — 1] as can be seen from their 
Fig. 1). Once the red and blue galaxy groups are defined 
(through the colour or SFR criterion), the authors compute 
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the mass functions of the blue and red groups and the con- 
tribution of blue and red galaxies to the total mass function. 
They show that the massive end of the mass function is dom- 
inated by red galaxies while the fainter end is dominated by 
blue galaxies. They then define two transition masses: 

(i) Mti is the mass for which the mass function of red 
galaxies equals that of the blue galaxies, i.e. both groups 
contribute equally to the total mass function. M tr shows 
only a moderate dependence on the red/blue separation cri- 
terion. It is found to increase from 10 10 5 M Q to 10 10 - 8 M o 
between 0.4 < z < 1.4. 

(ii) M quenc hB06 is the mass for which the fraction of blue 
galaxies drops below 1/3 (which means the contribution of 
blue galaxies is half the contribution of the red galaxies to 
the total mass function). The authors claim it represents the 
mass at which star formation "quenches" . It increases from 
1O 1O 73 M to 10 1123 M Q between 0.4 sC z < 1.4. 

Various other investigations consider different photomet- 
ric filters, colour cuts, SFR estimate s or mass estimates , 
but rely on simila r methodology (|Arnouts et al.l 120071 : 
iHopkins etaill2007l ). Given the complexity of the method 
and the number of steps and assumptions involved, it is 
clear that the physical meaning of the observed transition 
masses can be quite far from that of the transition masses 
we compute from our gas accretion measurements. However, 
it is admitted that they correspond to a real change in the 
evolution of SFR or specific SFR with galaxy mass. The de- 
tails are difficult to assess. For instance, does the transition 
mass signify a total shutdown of star formation or a smooth 
decline? Such details are bound to be fuzzy because we are 
looking at populations whose properties have an intrinsic 
dispersion, as they have different histories, environment etc. 
Having recalled these difficulties, we now proceed to make a 
number of remarks: 

(i) Using a stellar to dark matter mass ratio of 
M*/M DM = 50 we get M* ock w 1O 9 - 9 M at z = 2, which 
is significantly smaller than Mtt found in B06 at any epoch. 
Moreover, since M tr increases with z, the disagreement is 
bound to be stronger if we were to observe M tr (z = 2). A 
similar disagreement is seen for M str cam. 

(ii) On the other hand, we can define a mass M quenc h 
where the cold gas accretion rate on the central galaxy drops 
to zero, also corresponding to the mass where /cold — at 
the inner halo, instead of / co id = 0.5 as for M s hock' We find 
M qucnch = 1O 12 8 M , which translates to w 10 11 ,1 M Q in 
stellar mass. Although M que nchB06 is expected to increase 
between z — 1.4 and z = 2, the agreement is remarkable. It 
is also roughly consistent w ith the results of iPozzetti et al.l 
|2003t ); lFontana et al.l l|2004h . based on the K20 survey. This 
raises the question as to what we should adopt as threshold 
in /cold ■ Indeed we see that moving this threshold can result 
in a shift of more than a decade in the masses obtained. 
There is actually little reason for expecting / co id = 0.5 
haloes to match observed transition masses precisely. 

(iii) We find a difference between M s i 10 ck and M quen ch of 
about 1.2 — 1.3 dex. This is actually a measure of the sharp- 
ness of the transition from cold to hot accreting haloes. This 
contrasts with the finding of B06 where M tr — M quenc hB06 ~ 
0.2 — 0.4 dex, which suggests a much sharper transition. In 
future modelling works and observational studies, comput- 



ing the sharpness in mass of the transition from blue to red 
galaxies could help constraining the mechanism responsible 
for the transition. 

(iv) We find a quasi-constant M s h oc k, with only very 
slight evolution with z, compatible with the findings of ear- 
lier theoretical works, both numerical and analytical. In 
contrast, both Mtr and M quenc hB06 evolve strongly with z. 
On the other hand, it is difficult to check the evolution of 
M qU ench with redshift, since the only epoch where / co id ~ 
haloes exist in the simulation is z = 2. A naive linear extrap- 
olation of the evolution of / co id at larger redshifts suggests 
a quasi-constant M que nch, but this is simply impossible to 
check and would require simulating a larger box in order to 
get haloes massive enough to reach / co id = earlier in the 
life of the universe. 

(v) If we believe that M quen ch has a physical meaning 
similar to M quO nchB06 this apparent difference in their evo- 
lution with redshift is problematic. Moreover, a stellar to 
dark matter ratio M* /Mum evolving with z would make 
this issue even worse. Indeed, a constant M quenc h combined 
with a decreasing M /Mum with z would lead to a de- 
creasing M qU ench (z) , which would then be in even stronger 
disagreement with observations. 

(vi) On the other hand, M strea m does evolve with z, with 
a slope similar to that of M tr and M que nchB06- It would be 
interesting to compute the mass at which the cold fraction 
/cold drops to zero at the outer halo because such a mass 
may be comparable in absolute normalisation and in slope 
to M qU enchB06- However, this is difficult to even estimate 
from the current simulation, and the corresponding masses 
are clearly out of reach as can be seen on Fig. [4] 



Finally, it becomes clear that in order to compare our the- 
oretical transition masses to observed ones, one should use 
the same transition definitions in the simulation as in the 
observations. This involves building catalogues of galaxies 
(rather than dark matter haloes), computing their luminosi- 
ties and colours according to the age and metallicity distri- 
but ions of their stars u s ing st ell ar population models such 



Bruzual fc Chariot! (120031) : [Fioc fc Rocca-Volmerangel 
1 19971 ): iLe Borgne et al.l (|2004l ): ICoelho et al.l (|2007l ). The 



resulting colour-magnitude diagrams of the galaxy popula- 
tion could then be compared directly to observations suc h 
as the SDSS M r /(u - r) distribution jBaldrv et alj|2004h . 
This is the approach adopted in SAMs such as those of 
ICattaneo et alj (120061). The same colour cuts as in the obser- 
vations 1 Arnouts et al.ll2007l ; Ivan Dokkum et al . 2000) could 
then be applied and the corresponding red and blue galaxy 
mass functions constructed, along with the corresponding 
transition masses. Since the colour bimodality of galaxies 
falls into place only after z ^ 1.5 (the last redshift bin of 
lArnouts et alj (|2007l ) shows no bimodality in (NUV-r')/K 
at all), it should not be expected that the bimodality can 
be clearly seen in a colour-magnitude diagram even at the 
lowest redshift of the simulation. Moreover, since the AGN 
feedback now often proposed as the origin of the galaxy bi- 
modality has not been implemented in the MareNostrum 
simulation, it might never appear even if we could continue 
the simulation down to lower redshifts. In any case, the 
above colour-magnitude diagram synthesis and the corre- 
sponding colour cuts have to be carried out quantitatively 
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to check for the presence/absence of a galaxy bimodality in 
the MareNostrum simulation. 



6 CONCLUSIONS 

We used the Horizon-MareNostrum galaxy formation sim- 
ulation to study the processes involved in gas accretion on 
galaxies. We introduced mass accretion rate weighted statis- 
tics that allow us to quantify the mass accretion rates as a 
function of gas temperature and metallicity. 

Gas accretion is bimodal both in temperature and 
in metallicity, defining a hot and a cold accretion mode. 
The cold accretion mode is associated with a combination 
of metal poor filamentary accretion and dense metal-rich 
satellite galaxy disc surroundings, while the hot accretion 
mode features strong chemical heterogeneity, and a radius- 
dependent metallicity. 

We give an analytical fit to the metallicity of the hot 
accretion mode as a function of radius, which will be relevant 
for future SAMs. 

We define M s h oc k and M s tream, the halo masses for 
which cold and hot accretion contribute equally, at the inner 
halo and within the whole halo, respectively. Haloes more 
massive than M sn0 ck develop stable hot shocks, but may still 
possess cold gas filaments nourishing the galaxy disc. For 
halo masses larger than M s trcam, these filaments disappear. 
Afshock is found to be quasi-constant with z, while M s tream 
increases sharply. Our results for M s h oc k are in agreement 
with the analytical stability calculations of DB06, and our 
metallicity measurements support their original assumption 
at the inner halo. Conversely, our determinations of M s tream 
disagree with their original predictions. We show that as- 
suming a low metallicity for the filaments, as we measured 
(10 -3 Zq is similar to primordial abundances as far as cool- 
ing is concerned) brings their model in agreement with our 
results. 

This suggests that, as long as the stability of hot shocks 
and cold streams is assumed to be mainly driven by a com- 
petition between compression and radiative cooling, their 
analytical modelling is accurate within this model; it then 
depends critically on assumptions made about the metallic- 
ity, via the cooling function. 

We propose that in addition to the transition masses, 
two other observables should always be considered by future 
SAMs or numerical investigations of the origin of the galaxy 
bimodality: 

(i) the sharpness of the transition log(M qucnc h/M s hock), 
whose observable counterpart could for instance be the log 
ratio of the quenching mass to the transition mass as defined 
in B06. 

(ii) The evolution of transition/quenching masses with 
redshift. 

Modelling these two quantities should be helpful in pinpoint- 
ing the nature of star formation quenching in galaxies. 

Comparing the transition masses we obtain to ob- 
served transition masses is a difficult task, and we found 
only marginal agreement. The diffuse cold gas supply drops 
to zero at the inner halo for an estimated stellar mass 
M* cnch ~ fO 111 M0 at z = 2, which is remarkably close to 
the quenching masses observed by B06. Unfortunately, we 



are not able to constrain the evolution of M quen ch. However, 
we note that the evolution of the observed quenching mass is 
similar to the evolution of Mstream. In this respect, the agree- 
ment between measured and observed quenching masses sug- 
gests that simulating realistic galaxy populations does not 
necessarily require more ingredients than the physics already 
modelled in the MareNostrum simulations. We recall here 
that no AGN feedback has been taken into account in this 
work. To be more conclusive, better constraints on Mquench 
and its evolution are needed. These could be obtained by 
means of zoom-simulations of smaller boxes centred on mas- 
sive haloes to lower z. Since M quoncn is determined mostly by 
the few most massive haloes of the simulation, a larger box 
or more realizations of a box of the same size are required 
to improve the statistics. 

Conversely, additional ingredients may become neces- 
sary in order to prevent clumpy gas accretion from reach- 
ing the galaxy centre when it is significant (a study of the 
contribution of clumpy gas accretion will be carried out in 
a forthcoming paper). But they need not be in the form 
of AGN feedback. Although stripping of the hot halo of in- 
falling satellite galaxies is properly resolved in the MareNos- 
trum simulation, the interface instabilities cold disc/hot gas 
and cold filaments/hot gas are not. These could be the miss- 
ing physics. 

Similarly, some energy input into the hot component 
might be required to maint ain its temperature and avoid 
cooling flows at later stages. iDekel fc Birnboiml (2008) pro- 
posed that "gravitational quenching" could be a solution. It 
involves keeping the inner halo gas hot through interactions 
with cold dense gas clumps (drag), allowing to transfer the 
potential energy of these infalling clouds to the inner halo. A 
crucial step arising from our study is that of the interplay be- 
tween the hot gas bubble and cold streams/clouds. Are hot 
gas bubbles able to disrupt filaments connected to the disc 
via electronic conduction, turbulence, Kelvin-Helmhotz in- 
stability? Are cold streams immune to disruption thanks to 
increased pressure (feeling pressure of the hot gas) and thus 
higher density and consequently higher radiative cooling effi- 
ciency or the underlying dark matter stream? Can electronic 
conduction be impeded by local magnetic fields of suitable 
intensity? Is there a pressure/temperature/halo mass for 
which a given cold stream of given density/velocity/DM flux 
will be disrupted by these instabilities? It will be decades 
before cosmological simulations with the same size as the 
MareNostrum simulation will have the resolution required 
to resolve interface instabilities. Hence, as a first step, more 
restricted, idealised experiments are needed in order to in- 
vestigate these phenomena. 
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